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Abstract 

We outline the solution of a fundamental problem in quantum theory which has hitherto 
lacked a proper solution, namely finding the requisite quantum theoretical framework 
guaranteeing that the calculated inverse spontaneous emission rate of a moving atom, 
as a composite system of charged particles interacting with the Maxwell field, is slowed 
down exactly as in time dilation. 
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It has recently been shown that a neutral atom bearing an electric dipole moment 
moving in an external magnetic field can accumulate a quantum phase [1-5]. A mov- 
ing dipole may, under suitable conditions, exhibit a detectable Aharonov-Bohm phase 
shift [6] and the rotational motion of a Bose-Einstein condensate in a vortex state can 
induce a magnetic monopole [7,8] distribution or an electric charge distribution [8]. 
These and other effects associated with atomic motion continue to receive considerable 
attention and especially so with the advent of atom optics [9,10] and laser cooling and 
trapping [11]. At first sight it would appear that the requisite theory for the descrip- 
tion of pheneomena involving moving atoms could be constructed as a straightforward 
extension of non-relativistic quantum optics by incorporating the translational motion 
of the atomic centre of mass. 

In fact the need to incorporate the centre of mass motion in quantum optics theory 
had necessitated a re-appraisal of the corresponding non-relativistic quantum electrody- 
namical theory where investigations sought to elucidate how the division of the motion 
into centre of mass and internal motions is affected by the presence of the interaction 
with electromagnetic fields [12-15]. One of the main outcomes of these investigations 
was the emphasis on the role of the Rontgen interaction [16] energy term which cou- 
ples the electric dipole moment to an effective electric field involving the centre of mass 
velocity and the magnetic field of the light. 

It was Wilkens [17] who pointed out that a theory which excluded the Rontgen in- 
teraction would lead to spurious velocity-dependent effects when evaluating the sponta- 
neous decay rate of an excited electric dipole moving freely in electromagnetic vacuum. 
Wilkens extended his work to include the Rontgen interaction and evaluated the scat- 
tering rate into a given solid angle in a given direction, deducing that this was free 
of any spurious velocity-dependences [18]. He did not, however, proceed to ascertain 
whether the calculated total spontaneous emission rate based on his approach would 
be consistent with the requirements of special relativity. It turns out that Wilken's 
approach in fact leads to an incorrect result for the total rate. 



More recently, Barton and Calegoracos [19] highlighted the absence in the literature 
of a proper treatment of the quantum theory of spontaneous emission of atoms moving 
in a classically assigned trajectory. This is so even for the simplest case of a uniformly 
moving atom. They put forward a theory in which they considered a model system of 
an atom in which the nucleus is the centre of mass and merely provides the Coulomb 
potential binding it to the electron and this system was assumed to be interacting only 
with the scalar field as a simphfied representative of the Maxwell field. In addition to 
these simphfications, their theory called for a careful distinction between energies and 
Hamiltonians. As far as the authors are aware, a workable theoretical framework of the 
problem in which a real atom interacts with the full (vector) Maxwell field is hitherto 
unknown and it is our purpose here to furnish such a framework. 

Our theory makes use of unambiguous canonical techniques; Lorentz transforma- 
tions and a gauge transformation are two of its ingredients. It turns out that not 
only the Rontgen interaction features prominently in the theory, but that modified 
atomic internal energy levels and eigenf unctions, which produce an expected Doppler 
shift, play important roles in the characteristics of the emission process. Finally, one 
needs to distinguish carefully between projections of vector fields parallel and trans- 
verse to the atomic velocity vector. These different aspects of the problem conspire 
in a remarkable manner leading to the correct result, namely that the inverse total 
spontaneous rate of the moving atom follows the time dilation formula, as required by 
special relativity. 

The model of a real atom we consider here involves two oppositely charged particles 
of charges ei = —62 = e and finite masses mi and 1712. In the centre of mass frame 
(atomic frame) we denote the position vectors of the two particles by q'^ and and 
the electromagnetic scalar and vector potentials as 0'(r') and A'(r'), respectively. The 
Lagrangian density for the electromagnetic field including the interaction with the two 
charged particles is 



(A'(r') + V>'(r'))' - c\V x A'(r'))'l + J'(r').A'(r') - p'(r')0'(r') 




where the electric field is E'(r') = — (A'(r') + V'0'(r')), the magnetic field is B'(r') = 
V' X A'(r) and the current and charge densities of the particles are given, respectively, 
by 

J'(r') = e [q'Ar' - q',) - ^'^^(r' - q'^)] (2) 

p'(r') = e[5(r'-q;)-<5(r'-q',)] (3) 

The notation is such that the atomic frame (the rest-frame) is referred to as 5" and 
all quantities relative to this frame are primed. The laboratory frame is the unprimed 
frame and will be referred to as S, relative to which the atomic centre of mass moves 
at velocity R, and all quantities relative to S are unprimed. 

The Langrangian density in Eq.(|l|) can be recast in terms of the primed centre of 
mass coordinates, defined by 

(miq; + m2q^) _ , , 

^ - 1 q - q^ - q2 (4j 

where M = mi + m2 is the atomic mass. Note that in the primed frame, or rest-frame, 
we must have R' = 0, by definition. This allows us to carry out a Power-Zienau-WooUey 
gauge transformation [20] and straightforwardly obtain the new Langragian density 



[A'(r') + V<f)'{r') f - c\V' x A'(r'))'J - P'(r').(A'(r') + VV(r')) 

+M'{r').{V' X A'(r')) (5) 



where the polarisation and magnetisation vectors are expressed as full multipolar series 
in closed analytical forms 

V'ir') = E e J' ^A(q: - R') 5(r' - R' - A(q: - R')) (6) 

i=l,2 

M'{r') = E e J' dXX{4 - R') x q'j 6{r' - R' - A(q; - R')) (7) 

i=l,2 ^0 

The Lagrangian density in Eq.(|^) has a manifestly covariant form, viz 

£' = --F'^'T' - -G'^'T' (8) 



where in the primed S' frame F'^^ is the well known electromagnetic field 4-tensor [21] 
and G'^^ is the polarisation field 4-tensor [22]. Formally G'^'^ has the same form as F'^'^ 
but with the substitutions E' T*' and cB' — Al'/c. 

Lorentz invariance allows us to write the Lagrangian density in the unprimed (lab- 
oratory) frame S exactly as in Eq.(^, or Eq.(^), simply by removing the primes. The 
total Lagrangian in the unprimed frame can now be written by adding the familiar 
relativistic Lagrangian contributions from the two particles as follows 

L = -m,ql/-f{q,)-m2ql/i{q2)+ J d'r [| (^^(r) - c'B\r)) + P(r).E(r) + At(r).B(r) 

(9) 

where 7(g) = (1 — g^/c^)^^/^; the electric and magnetic fields are given by E(r) = 
— A(r) — V0(r) and B(r) = V x A(r). It is important to bear in mind that the 
unprimed polarisation and magnetisation fields 'P(r) and Ai{r) appearing in Eq.(H) 
are not those in Eqs.(|^) and (^). The primed polarisation and magnetisation fields are 
rest properties and the unprimed ones are related to them by relativistic connection 
rules involving a Lorentz transformation of the polarisation field 4-tensor G^i, [22]. The 
interaction Lagrangian density (the last two terms in Eq.(^) can thus be rewritten in 
terms of the primed polarisation and magnetisation by direct substitution as follows, 



int 



VUr') + 7 ( Pl(r') + X M'ir 



c2 



E(r) + 

A4j|(r') +7 (All(r') - R X P'(r'))] .B(r) (10) 

where the subscript || (±) denotes the vector projection parallel (perpendicular) to 
R and 7 = 7(-R). Note that in Eq. (p!oD the electric and magnetic fields have been 
left untransformed and we have emphasised the dependence on the space arguments r 
and r' which are connected by the Lorentz coordinate transformation equations. The 
term in Eg . (|T0|) involving the velocity, electric polorisation and the magnetic B field is 
identified as the Rontgen Lagrangian interaction term [12-16], while the term involving 
the velocity, magnetisation and the electric field is identified as the Aharonov-Casher 
term [8,23]. 



Having expressed the interaction Lagrangian in terms of the known rest atomic 
properties 7-*' and Ad' , we now turn to the particle Lagrangian terms (given by the 
first two terms in Eq. (^ ) and seek to express their sum in terms of the unprimed centre 
of mass and relative velocities, R and q, respectively, using the relations 

qi = R+^q; q2 = R-^q (n) 

Concentrating on the unprimed (S) frame we now make use of the fact that the internal 
dynamics of the atom are not affected by relativistic considerations other than through 
R (i.e. the motion of the electron round the nucleus is essentially non- relativistic). We 
may then expand the sum of the particle Lagrangian terms up to terms quadratic in q 
to obtain 

„2 



{miqf/l{qi) + m2qlh{q2)) ^ -MR" h + (q" + ^(Rq^' 



^ -Mc^ + \[mR' + ^^^(ql + ^^^)] (12) 
where /i = mim2/M is the reduced mass. 



After making use of Eqs.(|TO|) and (|T2|), the new Lagrangian emerging from Eq.(^ 
now becomes the starting point of the canonical procedure with R and q as the canon- 
ical variables for the atom and A(r) and 0(r) for the fields. The canonical momenta 
are P (conjugate to R), p (conjugate to q) and n(r), which is identified as — D(r), the 
electric displacement field, is the momentum conjugate to A(r), while the momentum 
conjugate to is zero. Since we are interested in spontaneous emission by an atom 
characterised by an electric dipole moment, we may now ignore magnetic interactions 
by setting all terms involving to zero. The final Hamiltonian emerging from the 
canonical procedure can be written as a sum of three terms as follows 

H = Hl + H'j + H,^t (13) 

where 



(15) 



int 



^ :.(P'(r') X B(r)) + (P'(r') x B(r)).; ^ 



(16) 



In Eq.(P^ is identified as the unperturbed atomic Hamiltonian and is seen to be 
clearly divisible into a centre of mass part and an internal part and we should note the 
appearance of the relativistic factor 7 in the latter. The potential U (q) in Eq. (|I4D is the 
inter-particle Coulomb potential in the unprimed (laboratory) frame. In the primed 
frame (rest-frame) the inter-particle Coulomb potential, denoted as U'{q'), arises in the 
multipolar formulation from an integral term containing the square of the irrotational 
part 7-*'^ of the polarisation field, together with infinite Coulomb self energies 

— / d^r' = U'(q') + infinite Coulomb self energies (17) 

On disregarding the infinite Coulomb self energies, one then transforms the inter- 
particle Coulomb energy U'{q') to obtain f/(g), the interaction in the unprimed frame. 
The simplest and most direct route is by following the force transformation argument 
[24] to obtain 

U'{q') 

^'^^ 7 47reo7g' ^ ^ 

The expression for ifj given in Eq.(^) is the familiar unperturbed field Hamilto- 
nian which can be quantised following the standard methods of quantisation for a free 
field in the laboratory frame. Finally, Hint-, given in Eq.(|16D is a new form of interaction 
Hamiltonian which couples the rest polarisation field to electromagnetic fields in the 
laboratory frame. Note the division of polarisation in terms of || and _L contributions 
and the appearance of the relativistic factor 7, explicitly and also implicitly via Lorentz 
transformation from r' to r. The last set of terms in Hint are identified as the Rontgen 
interaction, expressed in a symmetrised form. 



In order to describe the spontaneous emission process, we need to solve zero order 
eigenproblems for the fields and for the atom. The fields can be quantised by following 
the standard methods of quantisation based on the Hamiltonian in Eq.(|l^) in the 
laboratory frame. The solution of the atomic eigneproblem involves consideration of 
the Schrondinger equation if^\I^(R, q) = ii^\E'(R, q) which can be written as 



■2M^^ 27/i \dql ^ dql ^ ^i^dql I 47reo7g' 



vl>(R,q) = Evl>(R,q) (19) 



where, without loss of generality, we have taken the direction of the velocity R to be 
the z-direction. Equation (|T^) admits solutions of the form \I/(R, q) = e^^'^tp{q) where 
E = fi^K^/2M + e. Upon making the substitution q'^ = ^yq^ we obtain a Schodinger 
equation governing the internal states of a hydrogenic atom in the atomic frame (rest- 
frame) S' such that 76 = e' where e' are the internal eigenenergies in the rest-frame 
S'. 

The corresponding eigenfunction possesses the same formal expression in the two 
frames, but note that in the laboratory frame S the position vector would be q not q' 
which means that the atom will appear to Lorentz-contract in the direction parallel to 
the velocity, as should be expected. It is important to remember that these physically 
consistent modifications to the internal energy levels and eigenfunctions have only 
come to light because of the modified form of the internal kinetic energy term as 
well as the Coulomb potential energy terms appearing in Eqs.(|l4D and (piSf ). Without 
the asymmetry due to the presence of 7 in the internal kinetic energy terms and the 
dependence of the Coulomb potential energy on q', the familiar spatial symmetry of the 
hydrogenic Schondinger equation would have been lost, leading to angular dependence 
and, consequently, to spurious features arising from the lifting of the degeneracy of the 
energy levels. 

We are now in a position to consider the energy and momentum conservation ac- 
companying the process of spontaneous emission of a photon described in the unprimed 
(laboratory) frame as having wavevector k and frequency u when the internal energy 
of the atom changes from to ej. Conservation of momentum requires that we have 



Kj = Kj — k where Kj = Mli/ti is the initial centre of mass wavevector and Kj is the 
final wavevector in the laboratory frame. Conservation of energy, on the other hand, 
demands that we have 

where huj^ is the energy level difference in the primed (rest) frame S' and we have 
ignored the second order recoil energy. Note that Eq.(pO|) is equivalent to a Doppler 
shift in the photon frequency. 

Two cases in the calculation of the spontaneous emission rate will have to be 
considered relative to the laboratory (unprimed) frame S, namely (i) when the dipole 
moment vector is parallel to the velocity vector and (ii) when the dipole meoment 
vector is perpendicular to the velocity vector. If these two calculations yield exactly 
the same result, then the spontaneous emission is deemed to be isotropic i.e. free from 
angular dependence. Imposing the electric dipole approximation, P'(r') = d'5(r' — R') 
where d' = eq', we obtain for the transition matrix element squared, with Hint as given 
in Eq.(^) as the interaction, 

(21) 

where only transverse (i.e. solenoidal or divergence-free) electromagnetic fields are 
involved and we have written E instead of — Il/eo, anticipating free space quantisation 
[25] . The Rontgen term contains the average velocity of the atom before and after the 
transition due to the symmetrised term in Eq.(|l6D but we may ignore the momentum 
of the photon since this is small compared to the initial momentum of the atom. 

The free-space normalised electromagnetic fields can be obtained straightforwardly, 
remembering that we should identify two orthogonal wave polarisations. We choose 
the z-direction, i.e. along R, as the axis along which there will be either a magnetic 
field (i.e. transverse electric or TE) or an electric field (i.e. transverse magnetic or 
TM). We can then write for a given wavevector k = (/c, 9, 0) the following normalised 



\ME\ 



J'2 



Eii(R) 



+ 



7 



E,(R) 



2R 



X B(R) 



electric and magnetic field vector amplitude functions 

E™(r) = 7^ [cos(^) cos(0)x + cos(^) sin(0)y - sin(^)z] e'(^"-"*) 



Br(r) = (^)'^'[sin(0)x-cos(0)y]e^('^— 



Er(r) = (^) ' [sin(0)x - cos(0)y] e^(''--*) 



B^^(r) = — =^ [cos(^) cos(0)x + cos(^) sin(</.)y - sin(^)z] e'(^"-"*) 

(22) 

where is a normalisation volume. 

Turning finally to the Fermi golden rule formula in the unprimed (laboratory) frame 
S, we find that the spontaneous emission rate can be written as 

Ti ^ ^n^h ^ Jo Jo 'y^c'^{l - Rcos{9)/c)^ 

(23) 

where the sum over A denotes the summation over the two wave polarisations, TM and 
TE, and the right hand side emerges after performing the integration over k with the 
help of the delta function. Substituting for the matrix element squared from Eq.(pT]) 
and the electric and magnetic fields from Eq.(p2D, we find that the spontaneous emission 
rate for a dipole parallel and perpendicular to the velocity in the laboratory frame S 
are given, respectively, by 

" 47^(1 -i?x/c)4 7 

r = r' ^^ 3[(i + i?Vc^)(i + x^)-4i?x/c] ^ n .25) 

^ 873(1 - J?x/c)4 7 

where Fq = d''^ujQ^ /{S-Keofic^) is the free- space rate of spontaneous emission of the 

atom in the atomic rest-frame S'. Note that only the TM mode is involved in the 



evaluation of the parallel dipole rate and that both polarisations, TM and TE, are 
needed to obtain the perpendicular dipole rate. It is seen that there is no angular 
dependence, i.e. the rate of spontaneous emission is isotropic, and it does indeed vary 
hke a relativistic clock. 
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